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1. INTRODUCTION 

The baseline assembly complete configuration of Space Station Freedom is com- 
prised of an inner core and outer articulating parts as shown in Figure 1 . The attitude of the 
inner core, comprising the module cluster, center truss, etc., must be maintained close to 
a local vertical local horizontal (LVLH) attitude using active control devices such as Control 
Moment Gyros (CMGs) and/or Reaction Control System (RCS) thrusters. This attitude con- 
trol is necessary to accommodate microgravity and pointing requirements for user pay- 
loads located in the inner core. The power required for the Space Station is provided by 
photovoltaic (PV) arrays located on the port and starboard articulating structures. 

The PV arrays have to continuously track the sun to provide the required electrical 
power. Due to the motion of the Space Station along its orbit, variations in orbital geometry 
with respect to the sun, and the core attitude fluctuations with respect to the LVLH orienta- 
tion due to disturbing external torques, the orientation of the PV arrays must be constantly 
adjusted with respect to the inner core to track the sun. The physical rotary joints which 
perform this function are called the Solar Alpha Rotary Joints (alpha joints) and the Solar 
Beta Rotary Joints (beta joints). The alpha joints provide a relative rotational motion be- 
tween the inner core and the outboard truss. The beta joints perform the PV array orienta- 
tion adjustment with respect to the articulating port and starboard truss. Assuming the sta- 
tion is at an LVLH attitude, the alpha joint rotation rate would be the orbital rate. The beta 
joint rotation varies slowly over an orbit and the yearly variation follows the solar beta histo- 
ry. The definition of the solar beta angle is shown in Figure 2a. 

In this report, the authors investigate the power loss resulting from the elimination 
of one of the articulating joints. The objective of this study is to maximize available power 
with only one rotary joint. 

The motivation to eliminate the alpha joint is to reduce the cost of the station. The 
outboard truss and beta joint axis is aligned to minimize drag on the PV arrays and this cor- 
responds to a feathered array orientation for an LVLH core attitude as shown in Figure 1 . 
Optimal beta gimbal strategies are investigated. The power study is performed for LVLH 
and non-LVLH station attitudes. Intra-array shadowing, shadowing thresholds, and power 
loss from shadowing are also addressed. Alternatively, power loss due to elimination of 


i 


beta joints, with alpha gimbal reinstated is examined for comparison of available power. 


2. PROBLEM FORMULATION (BETA-TRACKING) 

Two optimal beta joint control problems are addressed in this paper and can be may 
be stated as follows: 

1) Evaluate the optimal variable beta angle history for the Space Station over a given orbit 
that maximizes the integral of the dot product of the unit vector from the spacecraft to the 
sun and the unit vector normal to the active side of the PV array over orbital daylight. The 
attitude of the core of the station is assumed to be LVLH. The dot product quantity is named 
direct sunlight factor in this study. The integral of the direct sunlight factor over orbital 
daylight provides a total measure of useful solar power over the entire orbit. 

2) Evaluate the optimal constant beta angle which has to be maintained over a given orbit, 
for the Space Station, that would maximize the same performance index as above. 

It is obvious that the solution to problem T ( pi ) would yield an equal or better per- 
formance index than the solution of problem 2 ( p2 ) since the optimal varying beta angle 
history could assume any value of beta. However, the solution to pi could introduce large 
rates for beta gimbal angles which may not be feasible. 

Once, the best beta gimbal strategy over an orbit is determined, the average power 
available over a year can be evaluated by incorporating the motion of the Earth around the 
Sun, and the nodal regression of the spacecraft orbit. The yearly average available power 
can also be computed for non-LVLH torque equilibrium attitudes (TEAs). Intra-array shad- 
owing and its resulting power loss can be obtained from PV array geometry and optimal 
beta tracking solutions and the results compared to shadowing losses for arrays with full 
sun-tracking. 


3. MATHEMATICAL MODEL 

The direct sunlight factor P can be written as a function of the orbital elements 
Q , / , ij, , y . , the 3-2-1 ordered Euler angles y> , 9 , 0 representing the station 

attitude from LVLH, the inner alpha gimbal angle a fixed at 90 degrees (as shown in Figure 
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1 ) to minimize drag and the beta gimbal angle control variable ft using coordinate trans- 
formation matrices. Here, Q represents the season dependent angle measured from the 
vernal equinox to the Earth-Sun line in the ecliptic plane and is related to the calendar day 
of the year by the expression Q = 2 n (DOY- 80)/365 . y is the inclination of the Earths 
equator to the ecliptic plane. </' is the orbit plane ascending node with respect to the vernal 
equinox and i denotes the orbit inclination with respect to the Earth equator, r] represents 
the true anomaly for the circular orbit. 

The coordinate systems (0,1,2) are shown in Figure 2 b [ref. 1 ] . The LVLH coordinate 
frame is as shown in Figure 2 a. The body coordinate system (4), the outer articulating truss 
coordinate system (5) and the PV array coordinate system (6) are as shown in Figure 1 . 
The alpha rotation is along positive y 4 or y 5 axis. The beta rotation is along positive z 5 or 
axis. The coordinate transformation matrices from Earth-Sun line coordinate system (0) to 
LVLH coordinate system (3) are as given below : 
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T* is the transformation matrix from LVLH to Body coordinate system (4) given by : 
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The transformation matrices from body to array coordinate system (6) are: 
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The unit outward normal vector from the active side of the PV array in array coordi- 
nate system (6) is [ 1 0 0| 7 , by definition. The Earth-Sun line unit vector [0 0 - 1] T can 

be represented in coordinate frame (6) by pre-multiplying it by transformation matrix . 
The direct sunlight factor P , which is the dot product of the above unit vectors can now 
be written as P = -7{j (1,3) during the orbital daylight and P = 0.0 during earth occupa- 
tion. Here, Tfj is a 3x3 transformation matrix defined by 7^j = 7lj * * rj * T\* T]* . 

In the initial portion of the study, it is assumed that the vehicle attitude is LVLH and hence 
the matrix 7^ is an identity matrix. 

Now the optimal control problem pi can be defined as follows: 

>h 

max [ P(fi) chf 

fi in) e ( rm I'/o- '//!) I 

4 



where, PWC denotes piece-wise continuous functions. One can chose tjq = 0.0 and 
ijf = 2 n without any loss of generality to represent one complete orbit. This is made equiv- 
alent to a typical Mayer type optimal control problem [ref. 2] using the the first order ordi- 

tlx 

nary scalar differential equation — = P(P) and initial condition r (/„) = 0.0 . 

chi 


The optimal control problem pi 


takes the form 

max x (ih) 
p (rj) E (PWC\t Mi \) 


subject to the differential equation — 


P(P) and initial condition x (/ () ) = 0.0 . 


Similarly, problem p2 can be written as 

max r (;/,) 

ft (tj) E (constant value P, f or all t] E [r/n, >//]) 
subject to the differential equation = P(fi) and initial condition r (fo) = 0.0 . 


4. SOLUTION PROCEDURES 

Since optimal control problem pi does not contain the state x on the right hand side 
of the differential equation, and since P is always positive semi-definite, the optimal control 
problem can be simplified to a parameter optimization problem where P has to be maxi- 
mized at every value of true anomaly during orbital daylight. The functional dependence of 
P on ft is obtained and the best value of p for a given // is obtained by satisfying the con- 
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For problem p2, a parameter optimization algorithm was used to obtain the optimal 
beta gimbal angle, which when kept constant over a given orbit maximized the perform- 
ance index. 

The optimal beta can thus be obtained for problem pi or problem p2 as explained 
above for a given orbit. For evaluation of performance over an year, a model of the regress- 
ing ascending node is used [ref. 3]. The daily rate of change of V 7 is as follows . 


dxfi _ 3 / 7(20 [ 

dt ~ 2(1 -e 2 ) a 2 


Cos{i) 


Here n is the orbital rate in radians/sec, and ('20 = -.001 0827 is the Earth oblateness coef- 
ficient. a v is the earth radius in meters while a is the orbit radius in meters, e denotes the 
eccentricity of the orbit and / indicates the orbit inclination as defined earlier. For a nominal 
altitude of 220 nm and inclination of 28.5 degrees, the rate of change of ascending node 
can be written as: 

^ _ 7.05078 degrees/day 

dt 

With this model and suitable initial values of the ascending node, yearly power losses due 
to optimal beta-tracking strategies can be investigated. 


5. NUMERICAL RESULTS (NO SHADOWING CON SIDERATIONS} 
a. LVLH Station Core Attitude : 


To clearly understand the variation in direct sunlight factor over an orbit, two extreme 
examples of orbits are used, namely, winter solstice (day 356) with ascending node 180 
degrees, and spring equinox (day 80) with ascending node zero as shown in Figure 2 c. 
The former example has a maximum magnitude of solar beta angle of 52 degrees while 
the latter example has a minimum magnitude of solar beta angle of zero degrees. 

Figure 3a depicts the variation of P over an orbit for winter solstice and ascending 
node 180 degrees. P has a value of unity for full sun-tracking during the orbital daylight. 


The region of occultation due to Earths shadowing is also depicted. Variation of P for opti- 
mal time varying beta angle (solution to pi), variation of P for optimal constant beta angle 
( solution to p2) and variation of the same with beta angle fixed at the solar beta angle are 
plotted. The time varying optimal beta angle solution is higher than the other two plots over 
the entire orbit as expected. The solar beta solution touches unity at at one point in the orbit, 
but has lesser average magnitude than the solution obtained by the optimal constant beta 
orientation. The percentage power available, defined as the ratio of the area under the plots 
of the single-axis beta gimbal strategies to the area under the full sun-tracking graph , was 
about 87.4% for optimal time-varying beta, 82.4% for optimal constant beta and 76.9% 
for beta gimbal angle equalling solar beta. 

Figure 3b similarly shows the direct sunlight factor P history for an orbit character- 
ized by a zero ascending node on the day of spring equinox. All the beta-tracking solutions 
are identically the same for most part of the orbit, except during the entry and exit from the 
Earths shadow, where optimal time-varying beta gimbal angle strategy produces a 1 80 
degree flip to obtain sunlight on the active side of the PV array. The percentage power avail- 
able was about 55.3% for optimal time-varying beta and 53% for optimal constant beta 
which was same as the solar beta. 

It was observed that the average power obtained for problem pi and the corre- 
sponding power for problem p2 were not very different for various orbits. The percentage 
of power improvement by using pi was in the order of 2-5% over using the constant beta 
gimbal angle strategy of p2. Figure 4a shows the beta histories for the orbit with ascending 
node of 1 80 degrees on winter solstice and Figure 4b depicts the same for the spring equi- 
nox, zero ascending node. The large rates (such as 180 degree flip in Figure 4b ) intro- 
duced by using pi solutions would not be feasible from an operational, attitude control, 
or structural dynamics point of view. Hence, in the forthcoming sections of the paper, only 
the optimal constant beta gimbal strategy is used. 

With the ascending nodal regression model as given earlier and with optimal con- 
stant beta gimbal tracking, a simulation was performed for one complete year with the initial 
ascending node assumed to be zero degrees at day zero. The ratio of the the available 
power using the optimal constant ( implies constant over an orbit, but slowly changes from 
orbit to orbit ) beta strategy to the power available using full sun-tracking is depicted in Fig- 
ure 5. The yearly average power available is about 60% of full sun-tracking power. The 
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yearly average power was re-computed for various initial conditions of the ascending 
node. However, the variation in average power was minimal and hence it can be tacitly as- 
sumed that the average available power is independent of the initial ascending node. The 
actual power available may be more due to reduced shunt losses when compared to full 
sun-tracking. The yearly variation of the optimal beta angle is compared to the solar beta 
angle in Figure 6, for an initial ascending node of zero degrees. The magnitude of the opti- 
mal beta angle is always greater or equal to the magnitude of the solar beta angle. This 
result has some special implication in intra-array shadowing, which is discussed in a later 
section. 


b. Non-LVLH Station Attitude : 

The steady state attitude of the core of the Space Station does not remain at LVLH, 
but oscillates about the Torque Equilibrium Attitude (TEA). The TEA is a function of the 
Space Station configuration, the operating altitudes and atmospheric disturbance torques. 
It is desired that the TEAs be within a ± 5 degrees about the LVLH, for payload pointing 
requirements. However, for earlier flight configurations, large pitch TEAs are possible [ref. 
4]. In this study, conservative maximum variations in the attitude of ± 10 degrees for roll 
and yaw and ±20 degrees for pitch are assumed. The only change in the analysis from 
the previous section is the value of the elements in matrix Tj . Yearly power availability (% 

of full sun-tracking power) is obtained as before for different TEAs in this range and the 
values are plotted in Figure 7. The variations in power availability are less than 2% for this 
range of attitude deviations from local vertical. 


6. NUMERICAL RESULTS WITH INTRA-ARRAY SHADOWING 

The analysis performed above assumes that the PV arrays do not shadow each oth- 
er at any time. However, the validity of this assumption depends upon the geometry of the 
arrays, namely, the width (2a) of the arrays and the distance ( b ) between the parallel beta 
axes and the solar geometry which is a function of the beta gimbal angle fi and the solar 
beta angle fi s as shown in Figure 8. The geometry shows the sun rays on the plane of the 
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paper which is a worst case shadowing phenomenon for the beta-only tracking problem 
( Figure 8 is always true for full sun-tracking arrays ). 

The model of the Space Station used, as shown in Figure 1 , has the inner PV arrays 
separated by 77 feet ( b = 23.47 m). The width of each array is 38.9 feet (2a = 1 1 .857 m). 
The port and starboard outer solar arrays are separated from the adjacent inner arrays by 
44feet(b = 13.41m). Figure 8 can be used to prove that for shadowing to occur the follow- 
ing equation must be satisfied; 


Sinfi 

tan(f -ft,) 


+ Cmfi > £ 

2a 


Since the inner array separations are larger than the separation between the adja- 
cent port or starboard arrays, the latter plays the dominant role in causing shadowing. Fig- 
ure 9 shows the variation of the left hand side of the above equation defined as the shadow- 
ing function with (i for various values of /?, . The right hand side of the above equation is 
a constant ( equal to 1 .131 17 for the two adjacent inner and outer PV arrays) for a given 
configuration geometry. Figure 9 clearly shows that for a given fi s , the shadowing function 
peaks at p = p s . Shadowing may occur only if the shadowing function curve crosses the 
b/2a line as shown in Figure 9. The smallest value of p s at which the shadowing function 

crosses this line for any ft is at fi x = 27.86 degrees. In other words, no shadowing occurs 
for any beta angle if the solar beta angle is less than 27.86 degrees. For solar beta angles 
greater than this value, the shadowing function cuts the b/2a line at two points. For a given 
solar beta angle, if the beta gimbal angle is such that the corresponding shadowing func- 
tion value is greater than b/2a , then partial shadowing between arrays will occur over some 
portion of the orbit. 

Figure 10 shows the typical variation in sun-tracking beta angle and constant opti- 
mal beta angle with the solar beta angle ( this is obtained from part of the data in Figure 
6 ) for an LVLH attitude. A shadowing threshold line is also plotted in the same figure. The 
shadow region is obtained from the beta gimbal angle at the intersection points of the shad- 
owing function with the b/2a line of Figure 9. The region of partial shadowing during por- 
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tions of each orbit is shaded in Figure 1 0. If the optimal beta angle strategy is used, shadow- 
ing occurs only for beta gimbal angles greater than the critical value of 35 degrees, which 
corresponds to the intersection of the optimal constant beta line with the shadowing thresh- 
old line. For full sun-tracking strategy, the beta gimbal angle is same as the solar beta (for 
LVLH attitudes) and hence shadowing will occur at angles greater than 27. 86 degrees. Prior 
to deployment of the outer PV arrays, no shadowing can ever occur, since the critical value 
corresponding to the larger b/2a is greater than the maximum possible solar beta angle 
of 52 degrees. 

The ratio of the the available power using the optimal constant beta strategy to the 
power available using full sun-tracking has been shown in Figure 5. A similar graph with 
shadowing considerations is shown in Figure 1 1 . The regions of occasional partial shad- 
owing using optimal beta tracking and the regions of continuous partial shadowing for full 
sun-tracking are shown. The former is obtained by shading the time regions in the beta- 
tracking plot, where the solar beta is greater than the critical value of 35 degrees. The latter 
is obtained by shading the time regions in the full sun-tracking plot, where the solar beta 
angle is greater than 27.86 degrees. 

The results indicate that, for the particular configuration studied, days of occasional 
shadowing constitute about 18% of the year for optimal beta tracking. Exact power loss 
due to shadowing can be obtained by more sophisticated ray-tracing methods, etc.. Days 
of continuous partial shadowing constitute about 30% of the year for full sun-tracking. 
Hence the actual ratio of available power for such configurations (which have shadowing 
problems) using optimal beta-tracking may be better than the nominal value of approxi- 
mately 60% as obtained in section 5 without shadowing considerations. 


7. A COMPARISON TO OPTIMAL ALPHA-ONLY TRACKING 

For alpha-tracking only, the alpha joints are free to rotate, and it is assumed that 
the beta angles are fixed at zero (as shown in Figure 1). The optimal alpha angle history 
to be determined, over a given orbit is that which maximizes the available power, assuming 
LVLH attitude. Occultation effects are also taken into consideration. The analysis is exactly 
the same as for problem pi , except the fact that the angle a (>/) is the control variable and 

fi = 0.0 is the constant beta gimbal angle. The zero beta angle is chosen since the solar 
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beta angle ( which is used for full sun -tracking beta gimbal angle ) varies from + 52 de- 
grees to -52 degrees. The Mayer problem as detailed earlier is constructed and the optimal 
control histories for alpha gimbal angle are obtained. 

Figure 12 shows the optimal alpha joint gimbal histories over one orbit for a best 
case of spring equinox with zero ascending node (zero solar beta angle) and a worst case 
of winter solstice with 1 80 degrees ascending node ( maximum solar beta angle of 52 de- 
grees ). The alpha history is linear, varying by 360 degrees per orbit. Figure 13 shows the 
variation of the direct sunlight factor over the two orbit cases described above. The direct 
sunlight factor is a constant over an orbit due to the fact that the normal to the orbit plane 
and the alpha gimbal axis are parallel. The direct sunlight factor is equal to the cosine of 
the solar beta angle corresponding to the particular orbit. Thus for the best case orbit, 
100% of the full sun-tracking power is obtained. For the worst case orbit, only 61 .57% of 
the full sun-tracking power is available. 

As performed in the beta-only tracking strategy, the model of the regression of 
ascending node is used to evaluate the average available power over one year. Figure 14 
gives the yearly variation of the ratio of power available from optimal alpha-only tracking 
to full sun-tracking power. This graph would be exactly the same as the variation of the co- 
sine of the solar beta angle over a year. The yearly average power available is approximately 
90% of the full sun-tracking power. This value is relatively independent of the initial ascend- 
ing node magnitude. 

Finally, the change in the yearly average power with variations in vehicle attitude is 
studied. Worst case plausible TEAs are chosen with variations in attitude of ± 10 degrees 
for roll and yaw and ±20 degrees for pitch. Figure 15 plots the available power (%) for 
these non-LVLH attitudes. The variations are minimal. 

It should be noted that the alpha-only tracking solutions with ft = 0.0 avoids any 
intra-array shadowing! If the Space Station configuration design has array geometry 
which causes shadowing for full sun-tracking, then the actual power loss by using alpha- 
only tracking may be even less than the average inferred value of 10%! 
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8. CONCLUSIONS 


a. Optimal Beta-tracking strategy : 

Assuming a 28.5 degrees inclination and LVLH attitude, optimal constant ( constant 
over an orbit) beta gimbal strategy produces an average of only 60% of the full sun-track- 
ing power over a year. On any given day the reduction may be as low as 17% or as high 
as 49% of the full sun-tracking power. Using an optimal time-varying beta gimbal angle 
improves the power availability by a small amount (about 2-5%) . However, this may require 
large beta angle rotation rates over an orbit. “Power reduction” has been used synony- 
mous with “sun-light reduction”. However, in reality, shunt losses are likely to be reduced 
for the optimal beta tracking when compared to full sun-tracking. 

For an LVLH attitude, an additional bonus for the optimal beta tracking strategy, is 
that the arrays are oriented edge-on (feathered) with respect to the velocity vector. This 
minimizes drag, increases orbital lifetime and decreases atomic oxygen degradation ef- 
fects. 

For variations in attitude of ±10 degrees in roll and yaw and/or ±20 degrees in 
pitch, less than 2% additional power loss occurs. 

For the configuration studied, prior to deployment of the outer solar arrays, no shad- 
owing between the inner port and starboard arrays were found. Subsequent to deployment 
of the outer solar arrays, occasional partial shadowing between adjacent arrays occurs for 
solar beta angles greater than 35 degrees. This constitutes about 65 days (18%) in an year. 
The shadowing only occurs for a few minutes in an orbit. For the same configuration, with 
full sun-tracking, continuous (through out the sun -lit portion of the orbit) partial shadowing 
occurs for solar beta angles greater than 27.86 degrees. This constitutes about 110 days 
(30%) in an year. 

b. Optimal Alpha-tracking strategy : 

Assuming a 28.5 degrees inclination and LVLH attitude, optimal alpha gimbal strate- 
gy produces an average of 90% of the full sun-tracking power over a year. On any given 
day the power availability may be as low as 62% or as high as 1 00% of the full sun-tracking 
power. 
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For variations in attitude of ±10 degrees in roll and yaw and/or ±20 degrees in 
pitch, up to an additional 1.5% power loss can occur. 

No power loss due to shadowing occurs for this strategy if the beta angle is kept 
constant at zero. 
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Figure 1 . Space Station Assembly Complete Configuration 


space station 



Figure 2 a : Definition of Solar Beta Angle 
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Figure 2 b : Definition of Intermediate Coordinate Systems 


Zero Ascending node orbit 



Figure 2 c : Definition of Best and Worst Case Orbits 



Winter Solstice (Ascending node = 180 deg) 
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Figure 3a : Variation of Direct Sunlight Factor over a Best Case Orbit 




Spring Equinox (Ascending node = 0.0 deg) 
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Figure 3b : Variation of Direct Sunlight Factor over a Worst Case Orbit 



Winter Solstice (Ascending node = 180 deg.) 
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Figure 4a : Beta Gimbal Histories over a Best Case Orbit 




Spring Equinox (Ascending node = 0 deg.) 



Figure 4b : Beta Gimbal Histories over a Worst Case Orbit 
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Figure 5 : Yearly Power Availability with Optimal 
Constant Beta Gimbal Rotation 



Initial ascending node = 0 deg 



Figure 6 : Yearly Variation of Optimal Constant Beta and Solar Beta 






Yaw angle (deg.) 
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Figure 7 : Yearly Power Availability with Optimal 
Beta Tracking for Non-LVLH Attitude 
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Figure 8 : Array Shadowing Geometry 
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Figure 10 : Beta Gimbal Regions of Solar Array Shadowing 
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Figure 1 1 : Yearly Power Availability with Regions of 
Occasional Shadowing 




Spring Equinox, Ascending node = 0 deg. (best case) 
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Figure 12 : Optimal Alpha History over Best and Worst Case Orbits 




Spring Equinox, Ascending node = 0 deg. (best case) 
Winter Solstice, Ascending node = 180 deg. (worst case) 
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Figure 13 : Variation of Direct Sunlight Factor over Best and 
Worst Case Orbits with Optimal Alpha Tracking 
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Figure 15 : Yearly Power Availability with Optimal 
Alpha Tracking for Non-LVLH Attitude 
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Maximum power availability while Beta-tracking is compared to full sun-tracking and 
optimal alpha-tracking . The results are quantified in orbital and yearly minimum, 
maximum and average values of power availability. 
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